Perturbative description of inclusive single hadron production at HERA 



S. Albino 

II. Institute for Theoretical Physics, University of Hamburg, 
Luruper Chaussee 149, 22761 Hamburg, Germany 

Abstract 

Light charged hadron production data in the current fragmentation region at 
^ - HERA are calculated using next-to-leading order perturbative calculations and 

Q ■ fragmentation functions obtained from similar data from e^e~ reactions. Gen- 

O ! eral good agreement is found at large photon virtuality and intermediate 

^ ' momentum fraction Xp, consistent with fragmentation function universality. 

The description of the small Xp and Q region is improved by incorporating 
, hadron mass effects. 

t> ; 

1 Introduction 

(~| , Unpolarized quark fragmentation functions (FFs) for charge-sign unidentified light charged hadrons 

Q-i| = vr^, and p/p have been well constrained by data from e+e^ + X using calcula- 

{2L|' tions at next-to-leading order (NLO) accuracy. Due to universality in the factorization theorem, such FFs 

^ . can be used to calculate the similar measurements of ep ^ e + + X. This contribution summarizes 

the main results of [IJ comparing ep reaction data in the current fragmentation region from the HI 
and ZEUS |3] collaborations at HERA with calculations using FFs extracted from e'^e~ reactions. 

'i> ' The kinematic degrees of freedom are chosen to be the centre-of-mass energy ^/s of the initial 

■ state ep system, the magnitude of the hard photon's virtuality = —q^, the Bjorken scaling variable 

Q ! X = Q"^ /{2P ■ q) and the scaled detected hadron momentum Xp = 2ph • q/q'^- The normalized cross 

^ I section (with the s dependence omitted for brevity) takes the form 

QQ . ^ 

O' , ,± /cuts '^^^'^^/x^'f rfr drrlO^ {x,Xp,Q'^) 

OO : FP™*- (cuts, XpA, XpB) = ^ r 7^2. dJ^^'^ ' 

I where "cuts" refers to a specified region in the (x,^^) plane, and where a;p^(5) is the lower (upper) 

^ ' edge of the Xp bin. The cross section and the kinematic variables are frame invariant, and aie measured 

. in the Breit frame, defined to be the frame where the photon energy vanishes. In this frame the target 

fragmentation region (xp < 0) contains the proton remnants, while the struck parton fragments into the 
current fragmentation region (xp > 0), and the latter process is equivalent to the fragmentation of a 
parton into an event hemisphere in e+e^ reactions. The factorization theorem dictates that, at leading 
twist, the highly virtual photon undergoes hard scattering with a parton in the proton moving in the 
same direction and carrying away an energy / momentum fraction y. One of the partons produced in 
this scattering undergoes fragmentation to the observed hadron moving in the same direction and 
carrying away an energy / momentum fraction z. In other words, after the change of integration variables 
z — > Xp/z and y x/y, the factorized cross section in the numerator of Eq. ([Hi takes the form 



daP;-t°^ _ dz dy ^ da'i ( x 

dxpdxdQ^ ^"^'"^^'"^ '~J.^ z J, y ^dzdydQ^ V^'"' ^^'"^^^ ^ 



(2) 

proton I ^ 2 \ j^h^ ( ,2^^ 



where fj"^°^°^ is the parton distribution function (PDF) of parton i in the proton, Dj"^ the FF of parton j 
to h^, da^^ the equivalent factorized partonic observable given to NLO in Ref. i^, iJ. the factorization / 
renormalization scale which distinguishes the soft from the hard subprocesses and as{fJ.'^) = ois{p) / (2vr). 



2 Comparisons with data 

At leading order in a^, eq. (dJ becomes 

^-»'- (cuts..,.,.,.) = j:,,iiQ^-)aAQ'^ . 

where the parton labels /, J are restricted to (anti)quarks qj only, which have electric charges e^j, and 
Gi{Q^) = J^^^^ dx xfj^"^"^{x, Q"^). In the limit that the Gi become independent of /, the numerator 
of Eq. Q is equal to the equivalent LO result for e+e~ — > + X, and therefore the two types of 
observables are distinguished only by the G/. If this discrepancy is small, a good description of HERA 
data is expected using EEs obtained from fits to e+e~ data, such as the AKK Q, Kretzer lH and KKP 
Q EE sets Q if universahty and fixed order (EO) perturbation theory are reliable. Calculations using these 
sets for the HI data are shown in Eig. [J] using the CTEQ6M PDE set |9| and the CYCLOPS program 




Fig. 1: Comparisons of theoretical predictions using tlie AKK, Kretzer and KKP FF sets witli the Xp distributions from HI l2l. 

IJO^ here and throughout this work, unless otherwise stated. The strong disagreement between the EE 
sets at large Xp most likely arises from large experimental errors on the EEs due to poor constraints from 
e+e~ reaction data at large momentum fraction. At high Q^, the calculation for all 3 EE sets agrees well 
with the data. Therefore, the disagreements at large and small Xp values found with the lower data 
may be due to effects beyond the EO approach at leading twist. Eor example, resummation of soft gluon 
emission logarithms that become large at small and large Xp may be necessary to improve the calculation 
here. This is illustrated in Eig. [2] (left) by the effect of scale variation on the calculation, being largest at 
small and large Xp. The effect of the observed hadron's mass nih is also important at small Xp for low 
values. Eor non-zero hadron mass, one has to distinguish between momentum, energy, light cone 
momentum etc., which are all equal when the hadron mass is negligible. According to the factorization 
theorem, the "momentum" fraction z appearing in eq. ^ is the fraction of light cone momentum carried 
away from the fragmenting parton by the observed hadron in a frame in which the spatial momenta of 
the virtual photon and the detected hadron are parallel, and Xp = ^p(l — fn\/ {Q'^^p)) should be replaced 
by the ratio of the hadron's to the virtual photon's light cone momentum, ^p. Using this approach [ ll, 
one finds that the experimentally measured quantity dcrP'^"*"'^ / dxpdxdQ'^ is related to the calculated 



'since this work was completed, 3 further sets |8| have been extracted using improved theoretical and experimental input. 



quantity daP™*""^ /d^pdxdQ^ by 



^^proto 



dxpdxdQ^ 



1 + 



^^proton h 

d^pdxdQ"^ 



(^5 Cp{Xp) , Q 



(4) 



which shows clearly that hadron mass effects become important at small Xp and low Q^. According to 
Fig. 121 this correction improves the description in this region, if we compare the results of this figure 
with the low results of Fig. [T] The choice rrih = 0.5 GeV represents an "average" mass for the 
light charged hadrons. We do not incorporate mass effects for the proton of the initial state, since this 
effect is expected to partially cancel between the numerator and denominator of eq. ([T]). By redoing the 
calculation with the MRST2001 PDF set lITTl . we see that the dependence on the choice of PDF set is 
small, particularly at small Xp, most likely because these quantities are well constrained but also because 
any variations in them are partially canceled between the numerator and denominator of eq. ([T]l- As for 
e~^e~ reactions, the dependence on the gluon FF is small, particularly at large Xp. 
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Fig. 2: As in Fig.[T] using only the AKK FF set. Left: The modifications arising from scale variation. Right: The modifications 
to the default predictions (solid line) arising from the replacement of the CTEQ6M PDF set by the MRST2001 PDF set of Ref. 
nil , from the removal of the evolved gluon, and from the incorporation of the hadron mass effect are shown. 

To further verify these observations and inferences, we perform similar calculations for the ZEUS 
data. The different FF sets lead to similar results and good agreement with the data at large and 
intermediate Xp (Fig. (3]). The scale variation (Fig. |4l top) generally decreases with increasing Q^, and 
is largest for small Xp. Both hadron mass effects and gluon fragmentation are most important at low 
and small Xp (Fig. HI bottom). 

3 Conclusions 

High measurements of inclusive single hadron production at HERA are well described by perturbative 
QCD in the framework of the factorization theorem using available FF sets. Although some disagreement 
is found with data at lower Q^, we note that there is significant room for improvement in the theory in 
this region, such as hadron mass effects studied in this work, but also resummation of the FO series at 
small and large Xp, higher twist effects and quark mass effects. Whether such effects are in fact relevant 
can be better verified by the effect of incorporating such data into global fits of EEs. More importantly, 
such data may also provide valuable information on the EEs' quark flavour components not constrained 
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Fig. 3: As in Fig.[T] but for tlie ZEUS data (3]. Eacli data set is measured in a specific a;p-bin and, togetlier witii its predictions, 
is siiifted upward relative to the one below by the indicated value for A. 



by e+e^ reaction data, particularly since these type of HERA measurements may now be made very 
precisely lfT2l . However, in order to constrain FFs for each hadron species individually, and to ensure 
that the data is not contaminated by light charged particles other than vr^, and p/p, the hadron 
species of the HERA data need to be identified. 
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Fig. 4: As in Fig.|2 but for the ZEUS data (H. 
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